A high redshift detection of the integrated Sachs— Wolfe effect 



O 
O 

(N 

CU- 
in 

(N 
> 
(N 
t> 
in 
r- 

o : 

VO . 
O ■ 

^ ' 

o ■ 
a: 



13 



Tommaso Giannantonio, 1 Robert G. Crittenden, 1 Robert C. Nichol, 1 Ryan Scranton, 2 Gordon T. Richards, 3 
Adam D. Myers, 4 Robert J. Brunner, 4 Alexander G. Gray, 5 Andrew J. Connolly, 2 and Donald P. Schneider 6 

'institute of Cosmology and Gravitation, University of Portsmouth, 
Mercantile House, Hampshire Terrace, Portsmouth P01 2EG, UK 
^Department of Physics and Astronomy, University of Pittsburgh, 3941 O'Hara St., Pittsburgh, PA 15260, USA 
3 Department of Physics and Astronomy, The Johns Hopkins University, 
3400 N Charles Street, Baltimore, MD 21218, USA 
^Department of Astronomy, University of Illinois at Urbana- Champaign, Urbana, IL 61820, USA 
5 Georgia Institute of Technology, Atlanta, GA 30332, USA 
''Department of Astronomy and Astrophysics, The Pennsylvania State University, University Park, PA 16802, USA 

We present evidence of a large angle correlation between the cosmic microwave background mea- 
sured by WMAP and a catalog of photometrically detected quasars from the SDSS. The observed 
cross correlation is 0.30 ± 0.14^7^ at zero lag, with a shape consistent with that expected for corre- 
lations arising from the integrated Sachs-Wolfe effect. The photometric redshifts of the quasars are 
centered at 2 ~ 1.5, making this the deepest survey in which such a correlation has been observed. 
Assuming this correlation is due to the IS W effect, this constitutes the earliest evidence yet for dark 
energy and it can be used to constrain exotic dark energy models. 



PACS numbers: 98.80.Es, 98.54. Aj, 98.70. Vc 



I. INTRODUCTION 

There is growing evidence that the expansion rate 
of the Universe is accelerating, which is believed to be 
the result of an unknown 'dark energy' which has come 
to dominate the present energy density. This is sup- 
ported by many measurements, particularly recent ob- 
servations of anisotropies in the cosmic microwave back- 
ground (CMB) anisotropies, such as by the Wlikinson 
Microwave Anisotropy Probe ( WMAP ) P, Q] , together 
with the present accelerated expansion inferred from the 
measurement of the Hubble diagram of the Type la su- 
pernovae H . Understanding the fundamental nature 
of this dark energy, what it might be and how it has 
evolved, is one of the biggest challenges facing cosmolo- 
gists. 

One way of probing the dark energy and its evolu- 
tion is through the integrated Sachs-Wolfe (ISW) ef- 
fect which produces anisotropies in the CMB rel- 
atively recently. While most of the CMB anisotropies 
were generated near the last scattering surface at the 
moment of recombination, when the universe was 400 ky 
old (z ~ 1100), additional anisotropies can be created 
later by gravitational interactions. If the gravitational 
potential <E> varies on large scales at some stage of the 
universe's evolution, the CMB photons will undergo an 
energy shift 
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where c = 1, r is the conformal time, the dot represents 
a conformal time derivative and the integration is along 
the line of sight of the photon. The ISW effect reflects the 
fact that a photon falling into a potential well will climb 
out at the same energy only if the gravitational poten- 
tial itself is constant in time, and will otherwise receive 



an energy shift depending on the evolution of the poten- 
tial. During the matter-dominated era, the gravitational 
potential remains constant and so $ = 0, leading to no 
ISW effect. However, if the Universe becomes dominated 
by curvature or dark energy, then new CMB anisotropies 
can be created. 

The observed features of the CMB anisotropy spec- 
trum indicate that the Universe is very close to flat, so 
a detection of the ISW can help constrain the dark en- 
ergy. Unfortunately, the ISW anisotropies will appear as 
a small addition to the anisotropies arising at higher red- 
shift; they cannot be easily distinguished because their 
signal is largest on very large scales where cosmic vari- 
ance is also large. One way to extract this signal is to cor- 
relate the entire CMB anisotropy map with some tracer of 
the dark matter distribution 0,00: the primary CMB 
anisotropies will not be correlated in any way with the 
matter overdensities we observe now, because all these 
structures have formed much later. 

The precise large scale CMB maps provided by the 
WMAP satellite have made observing such weak corre- 
lations possible. Many groups have detected the cor- 
relations between the WMAP CMB data and different 
tracers of the large scale structure, with results generally 
consistent with the predictions of the integrated Sachs- 
Wolfe effect for a dark energy model. The correlations 
have been detected with X-ray surveys 0, radio galaxy 
surveys H , infrared observations |l l| as well as op- 
tical surveys like the APM and the Sloan Digital Sky 
Survey (SDSS) [3 El El El 13 using both ordinary 
and luminous red galaxy samples. These surveys have 
also probed much different redshifts, allowing us to see 
the evolution of the correlation. 

In this paper, we use the NBC-KDE quasar catalog 
[H El from the SDSS. Objects in the NBC-KDE cat- 
alog are at z ~ 1.5, making this the highest redshift 
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sample ever used to probe the ISW effect. We find a pos- 
itive signal, suggesting that the dark energy behaves in 
a way compatible to the cosmological constant up to a 
redshift of 1.5. In the following we briefly describe the 
quasar catalog used (section[nJ , calculate its autocorrela- 
tion (section ITTTjl . Next we perform the cross-correlation 
with the CMB map (section lTV|) . Finally we describe the 
resulting constraints in section [V] followed by the con- 
clusions. 



II. THE QUASAR CATALOG 

The quasar data was derived from SDSS DR4 [H |U 

mmiSHIMIiilillilliil, usin S a nonparametric 
Bayes classifier method based on kernel density estima- 
tion (NBC-KDE) described in [HI- Briefly, this algo- 
rithm classifies quasars based on prior multi-color data 
on known quasars and stars, and is > 95% complete, 
with ~ 5% stellar contamination, to i = 21 [illlT& ]. The 
catalog contains N q = 344,431 objects with photometric 
redshift between 0.1 and 2.7 (see Fig. [!}, covering two 
distinct regions of the northern hemisphere of the Galaxy 
plus three narrow stripes in the southern, covering a total 
area of 6,670 square degrees. 

The stellar contamination is a potentially important 
systematic. Even if it does not contribute to the cross 
correlation (assuming the Galaxy has been cleaned from 
the CMB maps), it will still contribute to the quasar 
autocorrelation function. The stellar spatial overdensi- 
ties adds power on fairly large angular scales, which is 
difficult to explain with quasars alone. Even a stellar 
contamination as small as 5% can produce significant an- 
gular overdensities in the 5 — 10° range, where little con- 
tribution is expected from the quasars. For this reason, 
whenever modeling the expected behavior of our sample, 
we will assume that it is actually composed by a fraction 
k of stars and 1 — k of quasars. In section III we will 
show the contamination is k = 0.05 ± 0.01. 



Pixellation mask 

We are principally interested in the large scale corre- 
lations; to calculate these, we pixelize the quasar maps 
using the same HEALPix schema [3(J used to pixelize the 
WMAP maps. We perform most correlations with a res- 
olution parameter iV S ide = 64, corresponding to JV p i x = 
49,152 pixels of 0.92° resolution. Due to the partial sky 
coverage of the survey, only 16% of these pixels actually 
contain sources. 

Clearly the HEALPix pixelization will not exactly 
align with the SDSS regions where the quasars are ob- 
served. Some edge pixels will thus be only partially filled 
and it is important to take these effects into account: 
in the coarse pixelization described above, up to 20% of 
the pixels will be partially filled. To account for such 
effects, we use a high-resolution (iV s ^ s e h = 512) pixeliza- 
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FIG. 1: Redshift distribution of the quasars. A spline fit of 
this is used for the theoretical calculations. 



tion to determine the mask of the actual sky coverage of 
DR4, which will determine the fraction of each edge pixel 
that is matching the observed area. We base the mask 
on a random sample of galaxies in the DR4 database to 
ensure roughly uniform sampling in all directions. (A 
much larger stellar sample would be needed because of 
the high concentration of sources close to the Galactic 
plane.) By using a sufficiently large number of random 
galaxies (5 • 10 6 ) we can be sure to have good sampling 
when pixelized in the higher resolution. In the high- 
resolution map there are = 3 • 10 6 total pixels, of 
which only 5 ■ 10 5 cover the area of the survey, which 
means an average of 10 objects per pixel. In this way 
we estimate the coverage fraction of each low-resolution 
pixel, fi, as 



fi 
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(0 



64 



(2) 



where N^^ k (i) is the number of high-resolution pixels 
within the mask for each coarse pixel i, and for these 
resolutions, there are 64 high resolution pixels in each 
coarse pixel. 

We correct the maps by dividing the observed num- 
ber of quasars in a coarse pixel by the fraction of the 
sky within the pixel that was observed, yielding rij//j. 
For our correlation estimator, we down-weight such edge 
pixels by the fraction of sky they measure; this effec- 
tively accounts for the additional variance. A more con- 
servative approach is to simply drop these edge pixels, 
ignoring all quasars in them (although at this resolution 
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they contain roughly 20% of the catalog): we repeated 
our cross-correlation analysis using this schema and we 
found compatible results. 

We use the higher resolution to calculate the average 
number of quasars per coarse pixel, n: this quantity is 
the total number of quasars divided by the total area 
of pixels covered by the survey in the higher resolution, 
rcscalcd to the pixel surface area in the lower resolution: 



mask 

where N^^ k — ^maskM ^ s the total number of higher 
resolution pixels within the mask. 



Foregrounds 

There are a number of possible systematics in the cat- 
alog which could introduce errors resulting in a lack of 
completeness, bad redshift measurement or further stel- 
lar contamination; these could introduce artificial struc- 
tures in the maps and contaminate the measurements. 
We checked a number of these, including extinction by 
dust in our Galaxy, sky brightness, bright star obscura- 
tion and poor seeing in two different bands (r and g). 

The SDSS imaging data is obtained using drift- 
scanning, which produces long thin strips of data across 
the sky. Two adjacent strips are combined to make a 
stripe, which are then chopped into individual fields of 
dimension 10 x 10 arcmins [29|. Clearly, this observing 
strategy could introduce small correlations along a strip 
(or stripe), which could extend to very large angles (over 
100 degrees) in the imaging data, e.g. systematic dif- 
ferences in the zero-point calibration of the photometry 
in each strip. Such photometric calibration uncertain- 
ties were recently explored by [3l| and shown to be less 
than 2%, consistent over the whole SDSS area. This is 
below other statistical (shot noise) and systematic (ex- 
tinction, seeing) errors and therefore is not considered 
further here. We also note that the SDSS scanning strat- 
egy is not aligned with any cosmological or galactic signal 
and would therefore only introduce extra noise into our 
ISW detection rather than mimicing the signal. 

While the extinction is a quantity measurable for each 
observed object, sky brightness, seeing and number of 
point sources are global quantities of each 10 x 10 ar- 
cmins field of view. However, given an object we can 
find the foreground quantities associated to it through 
the ID number of its field of view and, because the fields 
are smaller (~ 1/25) than our pixels, we can consider the 
distribution of all these quantities in each pixel in the 
same way. 

From our random sample of SDSS galaxies, we find the 
value of each of these foreground quantities associated 
with each object (the extinction) or each field of view (sky 
brightness, seeing and number of point sources), and we 
build their distribution in each pixel. Then we take the 



median and we find the distributions of the medians of all 
pixels. Finally, we produce the masks for each foreground 
excluding the worst 20% pixels, i.e. the pixels whose 
median value for a given foreground is in the upper 20% 
tail of the distribution of the medians of that foreground. 
These masks are shown in Fig. [21 for the r band. 




FIG. 2: Foregrounds masks for extinction, seeing, point 
sources and sky brightness in the r band. The 20% of pixels 
with the worst contamination are shown in light green. The 
most relevant of these effects is the extinction of the sources. 



III. THE AUTO-CORRELATION 

To check the consistency of our method and to probe 
how biased our quasar sample relative to the underlying 
dark matter, it is useful to measure first the autocorre- 
lation function (ACF) of the quasar catalog. To do this, 
we use the estimator c**, where the index tt refers to the 
total catalog (including possible contaminations): 

where the sum runs over all the pixels with a given an- 
gular separation. As defined above, /j is the i-th pixel 
coverage fraction, n, is the number of sources in the i-th 
pixel, n is the expectation value for the number of ob- 
jects in the pixel. For each angular bin centered around 

N* = Y^fifi (5) 

is the number of pixels pairs separated by an angle within 
the bin, weighted with the coverage fractions. 

Here we present results using JVj, = 5 bins of in the 
range 0.5° < i? < 10°. We tried various angular binning 
schemes and the results seem fairly independent assum- 
ing a sufficient number of bins are used. Fig. [21 shows 
the ACF with and without the r band based foreground 
masks of Fig. [21 and the results are very similar using 
g band based masks. We find that the dominant effect 
is the extinction; the result obtained with this mask is 
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close to the one given by the application of all masks to- 
gether. This removal of the areas with the highest 20% 
of the extinction values is equivalent to cutting pixels 
with a reddening in the g band A g > 0.18, which is effec- 
tively what was done by |18| : for these reasons, we will 
use the reddening mask and not the others, in order not 
to excessively reduce the sample. We have also checked 
that a stricter cut in reddening (30%) does not change 
the result. For the other masks the 20% threshold is 
likely much more aggressive than required, but the in- 
dependence of the cross-correlation function (henceforth 
CCF) on these cuts shows that they are not significant 
contaminants. 



survey DR4 catalog; the quasar c qq (d) — b 2 c mm (d) is 
calculated from the matter power spectrum for the best 
fit WMAP third year model (WMAP 3), produced with 
cmbfast |35j | with a given source redshift distribution 
and assuming a linear bias factor, b relating the quasar 
clustering to the matter distribution. We have also to 
take in account the window function w($) associated with 
our pixelization, that is given by the HEALPIX team: 
the theoretical ACF c qq ($) is convolved with the window 
function «;(■#). The best values for the parameters are 
k = 0.05 ± 0.01 and b = 2.3 ± 0.2. The stellar contami- 
nation is thus in agreement with the expected value and 
both the stellar contamination and bias are consistent 
with those measured by 0, 0] . 
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FIG. 3: Auto-correlation function of the quasars measured 
for all the sample, for a single foreground mask and for all 
masks joint. A similar result is obtained for g band masks. 



This detection is consistent with the previous measure- 
ments [l8ll32l l33l l3~i| : these previous results used smaller 
data sets and were focused on smaller angular scales. In 
Fig. 0] we directly compare, using a similar binning, our 
detection to that of [13 which analyzed 80,000 objects 
from SDSS DR1 photometric catalog. 

We can model the total theoretical ACF c*' ($) as com- 
posed by the quasar and the star ACFs, c qq {d) and 
c ss (d), in the form 

c"(0) = (1 - k) 2 c qq ($) + fcV s (tf), (6) 

where k is the fraction of stellar contamination and we 
assume there is not any cross-term, due to the indepen- 
dence of stars and quasars (see ^H-) We obtain the stel- 
lar c ss (-d) from the average of 1000 subsamples of kN q 
stars (the number of stars we expect to have in the cata- 
log) from a random sample of 2 • 10 6 stars from the SDSS 



FIG. 4: The auto-correlation function of the quasar catalog 
with the reddening mask. The square (black) points are the 
observations c", the dashed line is the expectation c u , and the 
solid and pointed lines are its component (theoretical quasars 
and stellar contamination). We plot also the last points (red 
triangles) of the ACF measured by [Tg| for comparison. 



We have calculated the errors on the total ACF shown 
in Fig. 01 by producing 1000 random quasar maps, with 
the same statistics as the total catalog and an added 
Poisson noise. 



IV. THE CROSS-CORRELATION 

For the cross correlation analysis, we use the 
WMAP Internal Linear Combination (ILC) map derived 
from the third year WMAP data Q, pixelized in the 
same way and with the same resolution as the quasar 
map. Even though this ILC map was already built to 
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minimize the Galactic and other foreground contamina- 
tions, we have applied to it the most severe mask given 
by WMAP , the kpO mask, which corresponds to a cut 
of 32% of the sky. We checked that the results do not 
change significantly if we use the different frequency band 
maps V and W, corresponding respectively to 61 and 
94 GHz, all with the same masking; the results change 
slightly using the Q band map, which is the most affected 
by Galactic synchrotron contamination (see section [lV|l . 
We have also checked that the result remains consistent 
using the WMAP 1st year ILC map, and also does not 
depend on whether we use the smoothed or the raw single 
band maps. 

To measure the cross-correlation function (CCF) be- 
tween the quasar map and the WMAP ILC map, we used 
the estimator 

SAO -10 (£-«). <7) 

where Tj is the CMB temperature in the j-th pixel and T 
is the expectation value for the CMB temperature respec- 
tively. We again down-weight the partially filled pixels 
and Ntf is defined as above, but with a single weighting 
factor. We calculated this function in Nf, — 13 bins of d, 
in the range 0° < d < 12°, with and without using the 
foreground masks of Fig. [21 obtaining the results shown 
in Fig. |SJ We obtain very similar results using the r 
and g band masks. The reddening mask is the one that 
yields the lowest CCF; to be conservative and consistent 
with the ACF measure, we choose to apply this same 
mask. As expected, however, the reddening dependence 
is weaker for the cross-correlation measurement than for 
the quasar ACF. 

Fig. El displays the CCF between the WMAP3 ILC 
map and our NBC-KDE quasar sample. In reality, this is 
a measure of the cross-correlation between the CMB and 
a mixed sample of quasars and stars: although one does 
not expect a correlation between the cosmic radiation 
and local stars, we measured a small but non zero result. 
This indicates that the WMAP 3 ILC map, even after 
the most severe kpO masking, still has a small residual 
Galactic contamination. The stellar correlation has to be 
subtracted from the total detection yielding 

£ Ta (i>) . *"(«> - f W. ( 8) 
1 — k 

We compare this to the theoretical expected function 
c Tq (d) calculated again from the ACDM model with the 
WMAP 3 best fit parameters, using a program based on 
cmbfast [35|. 

To estimate the errors on the CCF and the covari- 
ance matrix, we use three different Monte Carlo methods. 
The first method is to produce a high number (2000) of 
random CMB maps with the WMAP best fit parame- 
ters and cross-correlate them with the true quasar map, 
after the application of the same kpO mask. Alterna- 
tively, we can do the reverse: use the true temperature 
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FIG. 5: Cross-correlation function of the quasars and the 
CMB measured for various foreground masks. Similar results 
are obtained for g band masks. 



map and make random maps of the quasars using the 
WMAP parameters and our observed bias, to which we 
add the Poisson error on the counts in the pixels. Both 
of these approaches produce similar answers, showing the 
covariances seem to be independent of any peculiarities 
of either of the two observed maps. 

These approaches give the covariances assuming the 
absence of correlations; while this should work well as- 
suming any true correlations are weak, it is important 
to understand the extent to which the presence of cor- 
relations will bias the covariance calculation. Indeed, 
if there are strong correlations, then these approaches 
should overestimate the errors. To account for corre- 
lations, we want to generate random temperature and 
quasar maps with the same ACF and CCF of the mea- 
sured maps, including also the Poisson uncertainty in the 
quasar counts. 

Based on the standard ACDM model, we can generate 
the expected angular power spectra for the anisotropics 
Cf T ,Cf q ,Cf q for the temperature only, the cross- 
correlation and the quasar autocorrelation. Here, the 
cross spectrum is assumed to arise solely from the ISW 
effect. From these power spectra, we can generate three 
random maps and use them to calculate the errors in the 
cross-correlation j3|| . We begin by making random tem- 
perature maps, T[, based on Cj T . (We neglect any noise 
which is thought to be small on the scales of interest.) 
We then decompose the quasar power spectrum into two 
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can now find the covariance matrix due to sample vari- 
ance, R^ mp . The results are consistent with what we 
obtain with the errors in the temperature only. 

The errors in the estimate made from Eq. © should 
also include measurement errors; assuming the mask is 
known, these can arise from uncertainties in k, the frac- 
tion of stellar contamination, or from the assumed stellar 
cross correlation c Ts . The full covariance is thus approx- 
imately: 

Rl^^ + tfRj'+wl (12) 

where 

J» = (cP - - £j s ) (13) 

and we have assumed the stellar contamination k <C 1 
and <7fc is its error. We account for the uncertainty of 
the star CCF Rjf calculating cross-correlations between 
random samples of kN q stars and the CMB map. The 
best fit and diagonal errors are shown in Fig. 



FIG. 6: The cross-correlation with the quasar catalog shows 
a small stellar contamination in the WMAP 3 ILC map with 
the kpO mask. The dashed line is the measured total CCF 
c Tt , while the point -dashed line is the measured stellar CCF, 
c Ta . The solid line is the difference between the two, which is 
our estimator for the true quasar CCF, c Tq . 



parts: 



Cf q = C," 11 + C, 



-,qq± 



(9) 



where the parallel and orthogonal signs indicate com- 
pletely correlated and uncorrelated with respect to the 
temperature map, and 
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(10) 



Using C 99 " and the same phases as for the temperature 

HI 

map, we create a correlated quasar density map, S i ; 
we add to this an uncorrelated quasar density map,^^ 
created using C? qJ ~ , with independent random phases. 



The total quasar density is i5[ 



§1 , and we can 



now build a random total quasar map n£, as 



1 



(11) 



Finally, we can add random Poisson noise to this, which 
we derive from the quasar number in each pixel. 

Generating 2000 Monte Carlo simulations n\ and cor- 
relating them with the random temperature map T[ we 
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FIG. 7: The measure of the cross-correlation. The points are 
the observed correlation between WMAP and the quasars, 
the solid line is the best fit ACDM theoretical model and the 
dashed line is the prediction for the WMAP 3 best fit model 
with b = 2.3. The points are highly correlated: the typical 
level of correlation between two neighbouring bins is ~ 95%. 



Three of the points in Fig. 0are more than lcr greater 
than zero, but the points are all highly correlated: the 
total significance of the ISW detection is found to be 
conservatively 2.1tr. Making less conservative choices in 
the analysis can lead to slightly higher significance, but 
it does not exceed 2.5a. This significance is based on 
a theoretical model for the expected ISW signal; using 
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a modified version of the cmbfast code, we calculated 
the predicted cross-correlation function for the WMAP 3 
best fit ACDM model, for a matter map with the redshift 
selection function shown in Fig. ^ We include the pixel 
window function and bin the expected correlation into 
the Nb — 13 Cj in the same way as we calculated the 
experimental CCF. 

We can compare the theoretical CCF Cj (the index Tq 
is understood) with the observed values Cj and assume a 
Gaussian likelihood model as 



C = (27r)- Ar / 2 [det%]- 1/ "exp[ 



1/2, 



E 

ij 



(14) 

where Rij is the CMB-quasar cross-correlation function 
covariance matrix defined in Eq. 1)12(1 . For a given dis- 
tribution of sources, the shape of the theoretical curves 
remains unchanged to a good approximation; however, 
the amplitude of the cross correlation strongly depends 
on the cosmological model, so that we can write [36| 



(ti)=A(Cl m ,w)g T «(ti) 



(15) 



where g T? ($) is normalized to 1 at 0°. We found the best 
value for A maximizing the likelihood, i.e. 



= l n ij 9i c j 
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(16) 



and the variance 
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FIG. 8: Frequency dependence of the cross-correlation, for 
the ILC map (black, solid), W (blue, long dashed), V (red, 
short dashed) and Q (green, dotted) bands. The first panel 
shows the measured CCF with the reddening corrected KDE 
catalog; in the second panel we plot the observed CCF with 
our random star samples (see text); the last panel shows the 
subtraction of the stellar contamination. 
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(17) 



The best fit for the CCF is A = (0.30 ± 0.14)/zK. 

In order to study the frequency dependence of our re- 
sult, we measure the cross-correlation using the single 
band CMB maps (W, V and Q); in the first panel of 
Fig. [8] we show the measured CCF for the single band 
maps, and we can see that the measure is almost fre- 
quency independent. We see in the second panel that 
the correlation with the random star sample is different 
for the different maps, being low as expected for the ILC 
map, increasing for the W and V bands, in which the 
Galactic contamination is more relevant, and being even 
bigger for the Q band, which is significantly affected by 
synchrotron radiation. The effect of the subtraction is 
shown in the last panel: the resulting CCF is still fre- 
quency independent and very consistent for the ILC and 
the V and W bands, while the galactic contamination 
starts to be important in the Q band, actually hindering 
the measure of the cross-correlation for this band. 



V. COSMOLOGICAL CONSTRAINTS 

To compare with these observations, we calculate the 
expected ISW cross-correlations based on linear theory 



using a modified version of the cmbfast code [37J. We 
calculate the quasar auto-correlations in the same way 
and assume a non-evolving linear bias factor. For the 
purposes of calculating the expected cross correlation, 
we use the actual measured redshift selection function 
of the sample ip = dN/dz normalized to unity. Were 
the quasar bias to evolve with redshift (e.g., as detected 
by EEl EES), this would effectively shift the red- 
shift weighting. We use the measurements of the quasar 
auto-correlation function to determine this bias to be 
b = 2.3 ± 0.2, consistent with previous measurements 
made at smaller scales [l8l l32l l33l 13^ . 

Assuming the cross-correlation that we see is due to 
the ISW effect, we can put some constraints on the na- 
ture of dark energy. First, consider the pure CDM model 
without any dark energy, (fl m — 1,Oa = 0); by relax- 
ing some assumptions, such as using a strongly broken 
power law for the primordial power spectrum, such mod- 
els might be consistent with the WMAP data [3^ . How- 
ever, these models would predict no ISW correlations, so 
would be disfavored at the 2a level with this data alone, 
and even more strongly when other ISW observations are 
included. 

Next, consider a flat dark energy dominated model 
(wCDM) with constant equation of state w. We first 
explore the likelihood function of the parameters fl m ,w 
with the constraint that the values of u>b = ^bh 2 ,Lo m = 
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£l m h 2 and the other parameters are fixed to the WMAP 3 
best fit values (ui b = 0.0223, w m = 0.128) [U- Here and 
below, the Hubble parameter is 100ft.kms~ 1 Mpc~ 1 . We 
obtain the result shown in Fig. |5J 
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FIG. 9: Likelihood contours within 1, 2 and 3 a on w — fl m 
inferred from this ISW detection alone. On the thick red line 
lies the family of the models whose TF spectrum matches the 
WMAP measured one, having the right comoving distance to 
the LSS, according to Eq. 1181 . Fhe thin black lines are the 
constraints on the Hubble parameter h {40| at I, 2 and 3 a 
(solid, dashed and dotted) assuming ui m = 0.128. 

We can see that the data are in favor of the 
ACDM model, but due to the weak detection only mod- 
els far away from this are actually ruled out. Most no- 
tably, models with a very small matter fraction predict 
too large a correlation and are inconsistent with the mea- 
surement. Fixing w — —1 yields the la interval for the 
matter fraction 0.075 < Q m < 0.475. 

While we have used WMAP constraints on the matter 
and baryon densities above, most models will actually 
be inconsistent with the positions of the CMB Dopplcr 
peaks. It is interesting to consider the family of mod- 
els which are consistent with the full temperature power 
spectrum measurements: for this, we need the angular 
scale of the Doppler features to be fixed to the observa- 
tions. The sound horizon scale is effectively fixed when 
we fix uj m and u>t>, so we must add the additional con- 
straint that the models have the same comoving distance 
to the last scattering surface, Df, given by 



#o Jx/ix+z.) v^ao 4 - 3 ^-) + a m a 

where z* is the redshift of the last scattering surface, 
weakly dependent on u> m . is kept constant if the 



variations in w are compensated by changes in the Hubble 
parameter h and the matter density f2 m : in Fig. [5] we 
show the family of the models fulfilling this condition, 
and we see that most of them are compatible with our 
ISW detection. 

This range of models is consistent both with the CMB 
autocorrelation and cross-correlation measurements: for 
instance, we show in Fig. HOI the temperature power spec- 
tra of two of these models; it is slightly different from the 
WMAP 3 best fit only at very large scales. 
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FIG. 10: Temperature power spectrum of models with the 
same angular diameter distance to the last scattering surface 
as the WMAP 3 best fit compared with the binned WMAP 3 
data (points). Changes in w are compensated by changes in 
fi m , but h undergoes considerable variations. 



Note however that many of these models are inconsis- 
tent with direct measurements of the Hubble constant. 
For w = —0.5, the Hubble constant would have to be as 
low as h = 0.55, while for w = —2.0 the Hubble constant 
would be unreasonably high, h = 1.20. Current limits on 
the Hubble constant, e.g. h — 72±8 |4jj, would constrain 
our measured w in the range —1.18 < w < —0.76. Mod- 
els with w in the range are practically indistinguishable 
from the best fit cosmological constant model plotted in 
Fig. El 

We have also investigated different classes of models 
which might be more likely to produce a significantly 
different ISW signal at the quasar redshift. There are 
two ways models might be strongly ruled out given our 
relatively weak detection: either they predict a correla- 
tion of the opposite sign, or they predict a much higher 
amplitude of correlation. Producing a negative correla- 
tion requires that the gravitational potential grow in time 
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rather than decay, which is difficult to arrange in typical 
dark energy scenarios because the accelerated expansion 
tends to slow down the growth of structure; one possibil- 
ity is a closed model without dark energy, as suggested 
by Nolta et al. [hJ. 

Producing a much larger signal at high rcdshifts is also 
difficult, given the other constraints on dark energy. For 
the models discussed above where dark energy scales as 
a power law, its fractional density tends to be small at 
a redshift of z ~ 1.5. In a cosmological constant model, 
^de(z — 1.5) — 0.16; while this density can be higher 
for w > — 1, the transition to dark energy domination 
becomes less sudden, leading to a smaller effect. If how- 
ever there were a sharp drop in the dark energy den- 
sity, it would be possible to be for the dark energy to be 
large at high redshifts while still remaining compatible 
with constraints from lower redshifts. In such a model, 
our measurement can limit the dark energy density at 
z = 1.5; models with Q,de{z — 1-5) < 0.5 would produce 
a much higher cross correlation than is observed, and can 
be ruled out at the 3a level. 

An alternative explanation of the dark energy prob- 
lem is to modify the laws of gravity on large scales; such 
theories may have consequences for structure formation 
which are significantly different to dark energy models, 
even if the background expansion appears the same. The 
ISW correlations are an important way of probing these 
differences, particularly at high redshift |4l|. There are 
many ways of implementing such changes, but much re- 
cent work has focused on the extra-dimensional model 
of Dvali, Gabadadze and Porrati (DGP) 142. |43l: under 
some assumptions, Song, Sawick and Hu |44| have re- 
cently shown the ISW signal to be comparable to dark en- 
ergy at low redshift, but significantly higher above z = 1 
in an open DGP model. At a redshift of z = 1.5, this en- 
hancement increases the expected cross correlation by a 
factor of two; while still consistent with our present obser- 
vations, larger quasar samples could be used to constrain 
such models in the future. 



Comparison of ISW detections 

We can now compare our detection with previous ones. 
Following previous convention [33, EH ; we plot the ob- 
served CCF at 6° in function of the mean redshift of the 
survey in Fig. ^2 This angular scale is chosen to avoid 
possible contamination from other effects that are domi- 
nant on smaller scales, such as lensing and SZ; however, 
it should be remembered that this representation is a one 
dimensional slice of the correlation function data. This 
approach can suppress the high redshift measurements, 
where a given angle corresponds to a larger physical scale. 
In the figure, we also add a conservative 30% error on the 
estimation of the mean redshift of the surveys. 

We also plot the theoretical expected values for the 
CCF at these redshifts for the models of Fig. ^| As 
above, we assume consistency with the CMB power spec- 



trum, i.e. we fix the comoving distance to the last scat- 
tering surface defined in Eq. I|18f) . while other parameters 
are fixed to the best fit WMAP 3 values. We see that 
the behavior is largely that expected from a cosmological 
constant model, with the amplitude dropping off at high 
redshifts. While many of the measurements are actually 
higher than expected, the differences are largely within 
the expected errors. This provides further support that 
the observed cross-correlations are due to the ISW effect. 
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0.1 1 
mean redshift (z) 

FIG. 11: Summary of the detections of the ISW effect 
through cross-correlation with different catalogs, compared 
with WMAP 3 best fit model. The blue (squared) points are 
in the order the correlations with 2MASS, APM, SDSS, SDSS 
high-z, NVSS+HEAO, as collected by |4^|: the green (trian- 
gular) points are the measure by |lq|. while the red (star) 
point is our KDE-QSO measure. The lines are the theoreti- 
cal expectations for WMAP 3 best fit model (solid) , and two 
models with w = —2 (long dashed) and w = —0.5 (short 
dashed) respectively (see text for details of the models). 



VI. CONCLUSIONS 

Here we have presented evidence of a weak correlation 
between the CMB and the distribution of high redshift 
quasars detected in the SDSS. Its amplitude, angular de- 
pendence and independence of the CMB frequency are 
all consistent with the interpretation as due to the inte- 
grated Sachs- Wolfe effect, with a significance in the range 
2 — 2.5<7, robust to changes in the mask and assumptions 
about stellar contamination. 

Without dark energy, no such correlation is expected. 
With a mean quasar redshift of z = 1.5, this represents 
the earliest evidence yet for dark energy and gives us a 
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means to further probe its evolution. Our measurements 
directly limit the density of dark energy at high redshifts, 
independent of its lower redshift behavior. They can also 
potentially provide interesting limits on alternative mod- 
els with modified gravity. 

These measurements will be improved when the pho- 
tometrically classified quasar data set is extended to the 
entire SDSS area. With a data set 40% larger, the photo- 
metric redshifts could be used to split the sample into two 
broad redshift bins, above and below z = 1.5, potentially 
allowing the evolution of the ISW effect to be seen within 
one self-consistent sample. To obtain even stronger cos- 
mological constraints, all the various ISW measurements 
should be combined carefully, including possible covari- 
ances which arise from the overlap of the different surveys 
in sky coverage and redshift; we are presently pursuing 
this |46j with an aim to provide an independent probe of 
the nature and evolution of dark energy. 
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